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Abstract. We present results from an individual particle 
based model for the collision, coagulation and fragmentation 
of heavy drops moving in a turbulent flow. Such a model 
framework can help to bridge the gap between the full hy- 
drodynamic simulation of two phase flows, which can usu- 
ally only study few particles and mean field based approaches 
for coagulation and fragmentation relying heavily on param- 
eterization and are for example unable to fully capture par- 
ticle inertia. We study the steady state that results from a 
balance between coagulation and fragmentation and the im- 
pact of particle properties and flow properties on this steady 
state. We compare two different fragmentation mechanisms, 
size-limiting fragmentation where particles fragment when 
exceeding a maximum size and shear fragmentation, where 
particles break up when local shear forces in the flow exceed 
the binding force of the particle. For size-limiting fragmen- 
tation the steady state is mainly influenced by the maximum 
stable particle size, while particle and flow properties only in- 
fluence the approach to the steady state. For shear fragmen- 
tation both the approach to the steady state and the steady 
state itself depend on the particle and flow parameters. There 
we find scaling relationships between the steady state and the 
particle and flow parameters that are determined by the sta- 
bility condition for fragmentation. 



1 Introduction 

Inertial particles in fluid flows have recently been subject 
of increasing in terest in sever al disci pl ines from dynam 



ical systems ("Benczik et al.l 
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by a flow. This problem displays already complex features, 
such as inhomo geneous spatial distributions (Maxey, 1 987t 
Wil kinson et al.l boovl) and multivalued velocity fields 



(IFalkovich et all I2OO2I: IWilkinson and MehhgL l2005b whose 



implications have yet to be understood completely. In most 
of these works a dilute regime is assumed, where particle 
collisions can be neglected. Some authors keep track of 
particle collisions numerically without considering the 
outco me of a collision, for exarnple to calculate collision 
rates dWang et al.[l200o(lBec et al.[l2005h . 



However, in many interesting applications, for example 
the growth of cloud droplets (Pruppacher and Klett, 1993) 
collisions of inertial particles play an important role. Pre- 
viously, this has mainly been studied using a mean field ap- 
proach, in the framework of which one treats the problem 
of particl e motion as a field eq uation. The Smoluchowski 
equation ( SmoluchowskiL [l917h is then used to model co- 
agulation and fragmentation of these particle concentration 
fields, instead of individual particles. Such an approach 
exhibits a number of problems. For example, the particle 
velocity may take on several values even at the same lo- 
cation of inertial particles when the dynamical attractor of 
the particle s fold s in the full velocity-position phase space 
fBec et aL , 2005 ) . Du e to the pre s ence of such 'caustics' 
(Fal kovichetalL I2OO2I; lBe3. 120031; IWilkinson and Mehhgl 
12005^ a mean field approach cannot be well founded and 
therefore relies on many assumptions and parameterizations. 
On the other hand, the simulation of multiphase flows us- 
ing full hydrodynamical models for each particle, e.g. with a 
finite element approximati on in an arbitrary Lagrange Eu- 
lerian framework (see e.g. iMaurvl Il999h can be very de- 
tailed but is usually restricted to an extremely small num- 
be r of particles (H igashitani et al., 2001; Zeidan et al., 200^). 
In IZahnow et al.l (i2008l) we therefore proposed a model for 
coagulation and fragmentation based on inertial particle dy- 
namics to help bridge the gap between the full hydrodynamic 
simulations and the me an field approaches. In the same year 
Wilkinson et al. IdZOOS") used a similar approach to model co- 



agulation and fragmentation of dust particles in an astrophys- 
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ical context. There smaU dust particles can gr ow into larger 
fractal clusters due to turbulent colUsions. In IZahnow et alJ 
(|2009a) the incorporation of such a fractal structure of the 
particle clusters in an inertial particle model was studied for 
marine aggregates. 

In this work we present results from such an inertial parti- 
cle based model for coagulation and fragmentation of heavy 
drops suspended in a synthetic turbulent flow as described 
in Section [2] Generally, we follow the model approach 
presented in IZahnow et al.l (12008. 2009b ) but focus here on 
the impact of particle properties and flow properties on the 
steady state size distribution of the drops that develops from 
the balance between coagulation and fragmentation. We 
study two different fragmentation mechanisms. First, par- 
ticles break up if their size exceeds a certain maximum al- 
lowed size. This is motivated by the hydrodynamical insta- 
bility o f hquid drops, for e xample rain drops settling due to 
gravity ( Villermauxl 2007 ). Second, particles fragment if the 
shear forces due to the fluid flow are too strong (see e.g. 

19991) . In contrast to our previous approach 



Thomas et al. 



we use here a lognormal distribution for the number of frag- 
ments, and compare this with other common approaches, 
such as binary or ternary fragmentation. 

Section [3] presents our results obtained from this model. 
We see that the distribution of particles as well as the mean 
average size in the steady state depends on the type of frag- 
mentation mechanism taking place. First, when fragmenta- 
tion occurs solely due to the particles exceeding a maximum 
stable size, the distribution is fairly uniform over all the ap- 
pearing coagulate sizes. Second, for fragmentation occurring 
under sufficiently large shear, the distributions typically de- 
cay exponentially beyond a certain coagulate size. In the case 
of shear fragmentation the mean average size of the coagu- 
lates depends strongly on the particle properties and the flow 
properties, such as the coagulate strength or the volume frac- 
tion of the particles. We show numerically that as a good 
approximation the influence of each parameter on the steady 
state can be treated separately and determine a decomposi- 
tion of the average coagulate size in the steady state with 
respect to the particle and flow properties. We find that the 
average coagulate size changes as a power law function of 
the particle and flow parameters, where the exponents of the 
power law are determined by the exponent in the stability 
condition for the shear fragmentation. 

The fluctuations over time in the coagulate size distribu- 
tion increase for larger mean average coagulate sizes. This is 
found to be a nonlinear effect resulting from the shear frag- 
mentation process. A similar decomposition and power law 
relationship with respect to the particle and flow properties 
as for the average coagulate size can be found here. 

For fragmentation due to particles exceeding a maximum 
allowed size the dependence on the particle properties is 
much weaker, instead the steady state depends mainly on the 
maximum stable coagulate size. However, the time to reach 
the steady state can vary greatly for different particle and 



flow properties. For example, flows with weak dissipation 
of turbulent energy allow for very large coagulates because 
there is almost no shear fragmentation, but since collisions 
also occur mainly due to shear the growth of coagulates is 
extremely slow. 

In Section |4] we show how in principle a scaling relation- 
ship for the average size of coagulates in the steady state can 
be derived, if size distributions, collision rates and fragmen- 
tation rates are known. Using approximate expressions for 
these quantities we again find a power-law relationship for 
the average coagulate size as a function of the particle and 
flow properties. Our calculations illustrate the dependency 
of the exponents of the power-law on the stability condition 
for the shear fragmentation. 

Section 12] gives a brief summary and states some conclu- 
sions. 



2 Mathematical Model 

In this section we present the mathematical model that will 
be the framework of our study. It contains a detailed de- 
scription of the dynamics of particles with inertia in a dilute 
suspension, a fairly general model for coagulation and frag- 
mentation of spherical drops. Additionally, the construction 
of a model flow is shown. We use a smooth random flow as a 
simple model for a turbulent flow below the dissipative scale. 

2.1 Dynamics of Dilute Suspensions 

In this paper we study suspensions of spherical inertial par- 
ticles of radius r, transported in an incompressible flow with 
dynamic viscosity /i. The particles are assumed to be much 
heavier than the surrounding fluid. We assume that the sus- 
pension is very dilute, i.e. particle-particle hydrodynamic 
interactions and feedback from the particles on the flow can 
be neglected. Additionally, we focus on a carrier flow with 
moderate Reynolds number and study only spatial scales be- 
low the Kolmogorov scale rj where the flow is sufficiently 
smooth. We therefore rescale space, time and velocity by the 
Kolmogorov length 77, time rj'^ /ly and velocity I'/rj (where v 
is the kinematic viscosity of the fluid). Assuming that the 
Reynolds number based on the particle size as well as the 
difference between the particle velocity V{t) and the flow 
velocity u{X, t) is small and the particle density is much 
higher than the density pp of the surrounding fluid, the mo- 
tion can be approximately described by the Stokes equation 
(.Max ey and Riley, 1983: .MichaehdeSi. 1 1 997i) . In dimension- 
less form it reads as 



(1) 



where X{t) is the d dimensional position of the particle and 
Strj — {2r^ ppv) / {^rf' ji) is the dimensionless Stokes num- 
ber. The effect of gravity has been neglected. 
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2.2 Coagulation and Fragmentation Model 

Next, we briefly describe the model for coagulation and frag- 
mentation that is us ed in this study. A mor e de tailed descrip- 
tion can be found in lZahnow et al. and lZahnow et al 
(l2009ab . 



The smallest particles considered will be called primary 
particles. These can combine upon collision to form larger 
particles, called coagulates. All coagulates are assumed to 
consist of an integer number a of these primary particles, 
i.e. the primary particles can never be broken up. A coag- 
ulate consisting of a primary particles has a radius r^") = 
Q,i/3j.(i)^ where 7-(^' is the radius of the primary particles. 
The coagulate's Stokes number depends on the radius, and 
tiierefore on a, with St^^^ = a^l^St^^^. Here S^^' is flie 
Stokes number for the primary particles. After the coagula- 
tion of two particles the velocity of the new particle follows 
from momentum conservation and the position is the center 
of gravity of the two old particles. To ensure that no colli- 
sions are missed, we use an efficient event-driven algorithm 
for particle laden flows (cf. Sigurgeirsson et all (I2OO1I) for 
details). Since hydrodynamic interactions between coagu- 
lates, that may affect the collision rates, are not included in 
such a model we approximate this by implementing a colli- 
sion efficiency Xc, which is the probability to coagulate upon 
collision. If particles do not coagulate upon collision they 
collide elastically. 

For fragmentation two different cases will be discussed, 
(i) Size-limiting fragmentation: If a particle becomes larger 
than some maximum number of primary particles Umax, it 
is broken up into k smaller fragments. For the ith fragment, 
where 1 < z < A: we set the new number of primary parti- 
cles ai to a random number drawn from a normal distribution 
centered around {olou ~ — i + 1) and with a 

standard deviation one, rounded to the nearest integer greater 
or equal to one. The last fragment contains the remaining pri- 
mary particles. This means that typically fragments will be 
of very similar sizes. This mechanism is motivated by the 
hydrodynamical instability of large water drops (e.g. cloud 
drops) settling due to gravity. We set the number of frag- 
ments to fc = 2 + ^, where ^ is a random number from a 
lognormal distribution with standard deviation one, rounded 
towards the nearest integer Such a distribution of fragments 
is a very common assumption for the fragmentation of drops, 
but later we will also comment on the implications of differ- 
ent choices for the number of fragments. 

(ii) Shear fragmentation takes place when the hydrody- 
namic force acting on the particle exceeds the forces hold- 
ing the coagulate together by a certain factor The hydro- 
dynamic force in this case is proportional to the shear force 



{2S,,S.,,tl^\ where 5., = \{M±M) 



S _ 
rate-of-strain tensor in the flow. iTavlor (11934 and later 
DelichatsiosI (Il975 ) derived an expression for the critical 
shear Scrit{r) across the drop required for fragmentation. 



under the condition that the characteristic time of drop de- 
formations is small compared to the time where this shear 
occurs. For our purposes this expression can be written as 



(1) 



(a) 



7a 



^1/3 



(2) 



where 7 is a constant, the coagulate strength parameter. If the 
shear force, calculated across the radius of the drop exceeds 
the threshold value given by Eq. (|2|i, the particle is broken 
up in the same way as for size-limiting fragmentation. 

The centers of the fragments are placed at a distance equal 
to the sum of their radii, perpendicular to the direction of the 
velocity and keeping the original center of gravity. The mag- 
nitude of the velocity remains the same to ensure momentum 
conservation. 

2.3 Model Flow 

We restrict ourselves to smooth, incompressible fluid veloc- 
ity fields since we focus on effects typically taking place on 
scales smaller than the Kolmogorov scale of a turbulent flow. 
To be able to perform long-term simulations at reasonable 
computational costs we consider synthetic turbulence in the 
form of a space-periodic, isotropic and homogeneous Gaus- 
sian random flow (Bee, 2005). Such flows are constructed to 
reproduce certain features of turbulent flows, but can not cap- 
ture all aspects of real turbulence, such as non-Gaussian tails 
of the velocity fluctuations and the energy cascade between 
scales. 

We write the flow as a Fourier series 



u{X,t) 



feez<'\{o} 



u{k, t)e 



(3) 



are the Fourier components, with the 
= u*{k,t) because u{X,t) is real- 



where u{k, t) e (i 
property u{~k,t) 
valued. The star denotes complex conjugation. By taking for 
u{k, t) the projection of a different vector v{k, t) e onto 
the plane perpendicular to the wave vector fc, incompress- 
ibility is ensured. The vector v{k,t) is assumed to be an 
Ornstein-Uhlenbeck process. It is a solution of the complex- 
valued stochastic differential equation 



dv ^ -^{k)vdt + (7{k)dW , 



(4) 



with ^(fe), <j{k) e M, where dW is a d dimensional com- 
plex Wiener increment. The parameters ^(fc), cr(fe) need to 
be chosen in such a way that the flow u{x,t) reproduces 
some features of a real turbulent flow, in this case the energy 
spectrum in the dissipative range of a turbulent flow. Here 
we use the exponential spectrum suggested by ICraichnan 

E{k) = C ■ {2TTkl/Lf exp{-P[2TTkl/L]) , (5) 

with j3 ~ 5.2 (see e.g. iMartinez et allll997l) and a suitably 
chosen normalization constant C. The constant I is the length 
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scale of coherent structures in the flow and L is the spatial pe- 



riod of the flow. We choose f (fc) — c and (T{k) = ^ cE{k). 
The constant c is then the inverse correlation time of the 
flow. The normalization constant is chosen in such a way 
that 2vk^E{k) sums to a desired value of the dissipation of 
turbulent kinetic energy e. The flow is then characterized by 
the correlation time 1 /c, the correlation length I and the dis- 
sipation e. 

If a fluid velocity field with few Fourier modes is cho- 
sen, no interpolation of the velocity at particle position is 
required, since it can be calculated from direct summation 
of the Fourier series. This allows for a resolution of the fine 




structures of the particle distribution in space. 



3 Coagulation and fragmentation in smooth random 
flows 

In this section we show results from the numerical simula- 
tion of the model described in section|2] we examine average 
quantities of the system and present size distributions for the 
different cases. 

3.1 Model parameters 

To reduce the computational effort we treat the case where 
the fluid flow depends only on two coordinates, i.e. we study 
a three dimensional flow where the velocity in the third direc- 
tion is negligible compared to the other two directions. We 
can then represent such a flow as two dimensional. Compar- 
isons were made with the full three dimensional case and no 
significant difference was found, except for a slowing down 
of the whole process due to a decreased number of collisions. 
For the fluid flow we take a total of 8 spatial Fourier modes 
in two dimensions into account, this is the lowest number 
for which isotropy is guaranteed. The period of the flow 
IS, L — 2tt. We set the correlation length of the flow to 
/ = 1, the correlation time to 1/c = 1 and the dissipation 
of turbulent kinetic energy to e = 1. This choice of para- 
meters results in large coherent structures in the flow (com- 
pare Fig. [U and a fast convergence to a steady state, due 
to sufficient collisions. The primary particles have a radius 
r(i)/L = 5 X 10-"^ and the Stokes parameter S*,'/^ = 0.05. 
For a 'typical' flow situation in a cloud this corresponds to 
primary particles in the range of 10~^m radius. We choose 
N = 10^ primary particles as initial condition. This implies 
a 2-d volume fraction of particles of approximately 0.08. 

The primary particles are uniformly distributed in the flow, 
with velocity V{0) = u{X{0),0). Due to the limited num- 
ber of primary particles, only a certain range of system para- 
meters is available. In particular the 'standard' values for the 
maximum stable coagulate size and the coagulate strength 
parameters are chosen such that in the steady state most pri- 
mary particles have formed larger coagulates, but the number 
of coagulates is still large enough to allow reasonable statis- 



Fig. 1. Snapshot of the position of 50000 particles with = 1. 
transported in a synthetic turbulent flow with dissipation e = 1. 
Due to dissipation the particles collect on a random attractor. 



tics. For size-limiting fragmentation the standard value is set 
at ctmax ~ 80 and for shear fragmentation the standard coag- 
ulate strength is 7 = 3.5. The standard value for the collision 
efficiency is Xc = 1- Unless mentioned otherwise, these are 
the parameters used to obtain the following results. 

3.2 Mean sizes 

We examine the numerical results for the full model with 
coagulation and fragmentation in terms of average quanti- 
ties such as the mean size of the coagulates. To characterize 
the steady state for the average number of primary particles 
per coagulate in the case of size-limiting fragmentation we 
could simply use the value for t ^ 00. This is not very 
precise, because the steady state is not static as shear forces 
fluctuate randomly over time. However, because the veloc- 
ity field is stationary, the mean number of primary particles 
in a coagulate {a{t)) will converge towards a constant value 

t+T 

ckoo = lim ^ Jds {a{s)) when averaged over a time inter- 

t-+oo J 

val T to remove random fluctuations. We use this quantity 
to characterize average coagulate sizes in steady state. We 
also compute the standard deviation CToo of the size distribu- 
tion as a measure of the width of the distribution. To remove 
random fluctuations (Too is also calculated as an average over 
a time interval T, in the same way as aoc ■ Here, we choose 
T = 100 for the averaging time, which was found to be a 
sufficiently long time interval to remove the fluctuations in 
the steady state. 

In the following we examine how these limiting values 
ttoo and (Too depend on the properties of the particles and the 
properties of the flow for the different fragmentation mech- 
anisms. The large parameter space makes it difficult to in- 
terprete results from the model. We therefore consider the 
sensitivity of the results to each of the parameters separately. 
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We restrict ourselves to the four most relevant parameters, 
namely the maximum stable coagulate size Umax (for size- 
limiting fragmentation), the coagulate strength 7 (for shear 
fragmentation), the collision efficiency Xc, the volume frac- 
tion of the particles, characterized by the total number N of 
primary particles and the dissipation of turbulent kinetic en- 
ergy e in the flow. 

3.2.1 Coagulate strength 

First, we examine the dependence of the average number of 
primary particles per coagulate on the maximum stable 
coagulate size a.max and the coagulate strength 7. These two 
parameters determine the binding strength of aggregates for 
the different fragmentation cases. 

Figure |2a) shows the results for the case of size-limiting 
fragmentation, and a^o both increase with the maximum 
stable coagulate size Umax- Here, we find that aoo oc amax 
and (Too Umax- The proportionality constant is determined 
by the distribution of fragments during breaking, a fit gives 
a proportionality constant of approximately 0.5 for ctoo and 

0. 25 for (Too. The proportionality constant depends on the 
details of the fragmentation mechanism. For example, the 
case of binary fragmentation, i.e. splitting into two fragments 
shows the same scaling but proportionality constants of ap- 
proximately 0.66 and 0.2, respectively. 

Figure |2jb) shows the results for the case of shear frag- 
mentation. A first qualitative estimate of the shape of this 
Oicritil) curve can be derived from a scaling argument, 
which was already mentioned in iZahnow et alJ (l2008h . Solv- 
ing Eq. |2]for a given value of the shear in the fluid we obtain 
a critical coagulate size for this shear This critical coagu- 
late size is proportional to 7-^. We therefore expect that aoo 
scales the same way. This is indeed close to the result of 
the numerical simulations, where we determine a relation- 
ship 000(7) oc 72-6±o i (dashed line in Fig. I2b)). A more 
detailed theoretical argument for this scaling will be given in 
Sec.g] 

However, since the shear in the flow fluctuates over space 
and time there is no single critical size for coagulates. We 
therefore expect that the width of the size distribution will 
depend, among other factors, on the fluctuations of the shear 
in the flow. From Eq. [2] it follows that larger coagulates are 
more sensitive to fluctuations in the shear This can be seen 
by considering how a change of the shear from S" to 5* + AS* 
changes the value of Ucrit- We obtain 

Aae„t = 7' {{S + A5)-3 - {S)-^) , (6) 

1. e. the fluctuations in the value of acrit are expected to 
increase proportionally to 7^^. We therefore expect that the 
width of the size distribution will also increase proportional 
to 7'^. This is again similar to the result of the simulations, 
where we find (Too (7) oc 72-6±o.i (dotted line in Fig.|2b)). 
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Fig. 2. Variation of the binding strength of coagulates, (a) For the 
case of size-limiting fragmentation the average number of primary 
particles per coagulate (triangles) and width of the size distribution 
in the steady state (squares) as a function of the maximum stable 
coagulate size amax- The fits are 0.5amax for the average (solid 
line) and 0.25amax for the standard deviation (dashed line), (b) 
For the case of shear fragmentation the average number of primary 
particles per coagulate (triangles) and width of the size distribution 
in the steady state (squares) as a function of the coagulate strength 
7. The fits are 2.27^ ^^ for the average (solid line) and 2.07^ ^^ for 
the standard deviation (dashed line). All error bars are obtained 
from an ensemble of 10 different simulation runs. The inset illus- 
trates the scaling behavior of the average and the width of the size 
distribution, the dashed line is the mean. 



3.2.2 Collision efficiency 

Second, we examine the influence of the collision efficiency 
Xc on the average number of primary particles per coagulate. 

For size-limiting fragmentation the simulations indicate 
(Fig- Etb)) that in this case there is almost no dependence 
of the average number of primary particles per coagulate in 
steady state on the collision efficiency. Both aoo and the 
width of the size distribution a^c remain almost constant with 
varying Xc- However, it should be noted that while the col- 
lision efficiency does not seem to have a large impact on the 
steady state, the transient behavior is greatly influenced by 
the collision efficiency. In particular, the time to reach the 
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Fig. 3. Average number of primary particles per coagulate (trian- 
gles) and width of the size distribution (squares) in the steady state 
as a function of the collision efficiency Xc i e. the probability to co- 
agulate upon collision in the case of (a) size-limiting fragmentation 
and (b) shear fragmentation. The fits are 54.62Xc'^^ for the average 
(solid line) and 44.7xc'"^^ for the standard deviation (dashed line). 
All error bars are obtained from an ensemble of 10 different simu- 
lation runs. The inset illustrates the scaling behavior of the average 
and the width of the size distribution, the dashed line is the mean. 



Fig. 4. Average number of primary particles per coagulate (trian- 
gles) and width of the size distribution (squares) in the steady state 
as a function of the total number of primary particles A'^ in the flow 
in the case of (a) size-limiting fragmentation and (b) shear fragmen- 
tation. The fits are 1.37V°'^ for the average (solid line) and l.liV'''^ 
for the standard deviation (dashed line). All error bars are obtained 
from an ensemble of 10 different simulation runs. The inset illus- 
trates the scaling behavior of the average and the width of the size 
distribution, the dashed line is the mean. 



steady state increases greatly with decreasing Xc, for both 
size-limiting and shear fragmentation. 

For shear fragmentation an increase in collision efficiency 
increases both the average number of primary particles per 
coagulate ctoo and the width of the size distribution doo (Fig. 
Ob)). The numerical results suggest a dependency of the 

form Ofoo, CToo OC ;^;0.31±0.03^ 

This increase of the average and width of the size distri- 
bution with increasing collision efficiency can be understood 
in terms of the balance between coagulation and fragmenta- 
tion. In the steady state the size distribution, and therefore all 
its moments, including the average and the width of the size 
distribution are determined by the balance between coagu- 
lation and fragmentation. Increased coagulation due to in- 
creased collision efficiency requires a corresponding increase 
in the fragmentation, which in turn requires larger coagu- 
lates. Unfortunately, deriving an equation for these moments 
of the size distribution from this balance condition requires 



the apriori knowledge of the shape of the size distribution 
and equations for the collision and fragmentation rates. It is 
therefore not a trivial task. In section]?] we show how such a 
calculation can be carried out if size distributions, as well as 
collision rates and fragmentation rates are known. This cal- 
culation will formalize the above argument and show more 
clearly how the scaling of the steady state with the particle 
and flow properties can be understood. 

The reason why such a scaling does not happen in the case 
of size-limiting fragmentation is that this specific fragmenta- 
tion rule serves as a 'brick wall' for the size distribution. For 
all coagulates below the critical size the fragmentation proba- 
bility is zero, for all coagulates above it it is one. Therefore, a 
balance between coagulation and fragmentation is only pos- 
sible for one specific size, independent of how coagulation is 
increased or decreased. 
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3.2.3 Volume fraction 

The third important parameter of the system is the volume 
fraction of particles in the flow. This is defined as the ratio 
of particle volume to fluid volume. The volume fraction is a 
function of both the total number of primary particles in the 
flow N and of the radius ri of the primary particles. Since 
the impact of varying ri is very similar to that for varying TV 
we only focus on the variation of the total number of primary 
particles (FigurelDl. 

For size-limiting fragmentation we again find almost no 
dependence of the steady state on the initial number of pri- 
mary particles (see Fig. Ufa)), only the time to reach the 
steady state decreases with increasing number of primary 
particles. For shear fragmentation our results show that 
the average number of primary particle per coagulate in the 
steady state as well as the width of the coagulate size distri- 
bution increase with N (see Fig. Sfb)). The numerical results 
suggest a relationship of the form aoojCoo oc 7V° '^*° °'^. 
This dependence on the number N of primary particles in the 
system can be understood in the same way as for the collision 
efficiency Xc, since both parameters increase the coagulation 
probability in the system. We will also illustrate this in the 
calculations in Section]?] 

3.2.4 Dissipation of turbulent kinetic energy 

Finally, we examine the impact of the flow on the simulation 
results. The mixing properties of the flow can be easily var- 
ied by adjusting the dissipation of turbulent kinetic energy e 
in the flow. A higher e results in higher shear in the flow and 
therefore influences both the coagulation and the shear frag- 
mentation. Coagulation increases with increasing e, since the 
collision rate due to shear increases but shear fragmentation 
also increases with e. 

For size-limiting fragmentation, again we find almost no 
dependence of aoo and aoo on the parameter (Fig. ]5ja)). The 
steady state seems to be independent of the changes in the 
flow. For both fragmentation rules the transient behavior is 
strongly influenced by the value of e. For small e collision 
rates in the flow are very low and the system takes a very 
long time to reach the steady state. For example, for shear 
fragmentation with e = 0.1, the system is still in a transient 
behavior for t — 500, while for e = 1 the system reaches 
its steady state at approximately t = 25. In the case of 
shear fragmentation, the average number of primary particles 
per coagulate in the steady state is proportional to g-i-2±o.i^ 
and the width of the size distribution also decreases approxi- 
mately proportional to g-i-2±o.i (pig ]3b)). This means that 
fragmentation increases faster with e than coagulation. 

3.2.5 Discussion 

In summary, we find that for size-limiting fragmentation the 
average number of primary particles per coagulate in the 
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Fig. 5. Average number of primary particles per coagulate (trian- 
gles) and width of the size distribution (squares) in the steady state 
as a function of the dissipation of turbulent kinetic energy e in the 
flow in the case of (a) size-limiting fragmentation and (b) shear frag- 
mentation. The fits are 52.92e~^'^ for the average (solid line) and 
29.04e~^ '^ for the standard deviation (dashed line). All error bars 
are obtained from an ensemble of 10 different simulation runs. The 
inset illustrates the scaling behavior of the average and the width of 
the size distribution, the dashed line is the mean. Stars denote the 
respective quantities minus their offset. 

Steady state and the width of the coagulate size distribution 
are mainly determined by the maximum stable coagulate size 
Oimax- For the fragmentation mechanism used here we find 

aoo ~ 0-5amax and (Too ~ 0.25a,„Qa;. 

For shear fragmentation the average number of primary 
particles per coagulate in the steady state and the width of 
the coagulate size distribution vary greatly with all system 
parameters. We obtain a relationship of the form 

aoo,CToo (X7^'Xc'^^'e^S (7) 

where the scaling exponents Xi are given byAi ~ 2.6±0.1, 
As = 0.31 ± 0.03, A3 = 0.3 ± 0.03, A4 = -1.2 ± 0.1. 

Further simulations confirmed that the scaling relationship 
for each parameter can indeed be approximately determined 
independently of the value of the other parameters. This 
means that Eq. (]7]i is expected to be valid for this model for 
all reasonable values of the system parameters, i.e. parame- 
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Fig. 6. Scaling exponents \i for the different parameters 
{'y,Xc,N, e) as a function of the inverse exponent for the critical 
shear 1/^ (Eq. l[2j). The scaling exponents are linear functions of 
1/^. (a) Scaling exponents Ai for the coagulate strength 7 and A4 
for the dissipation of turbulent kinetic energy in the flow e. The fits 
areAi = 0.273 + 0.763- i and A4 = 0.091-0.429- i. (b) Scaling 
exponents A2 for the collision efficiency Xc and A3 for the number 
of primary particles A'^. The fits are A2 = 0.126 + 0.061 ■ \ and 
A3 = 0.117 + 0.063 - i. 



ter values that lead to 1 ^ aoo <C N . We note that both for 
increasing N and x the average coagulate size scales with an 
exponent of 0.3 ± 0.03. This indicates that it is equivalent to 
vary the number of particles or the collision efficiency since 
both influence the average coagulate size in the same way. 

When looking at the exponents A^ for each parameter, one 
can ask the question whether these are related to the specific 
form of the stability condition Eq. (|2|i and in particular to the 
exponent appearing in this equation, as the simplified argu- 
ment for the dependency of aoo on 7 in Sec. |3.2.1| suggests. 

To examine this connection between the scaling exponents 
Ai and the exponent of the fragmentation condition in Eq. 
(|2]l, simulations with different exponents in Eq. (|2]l were per- 
formed. For example, changing the exponent for the criti- 
cal shear to —1/2 instead of —1/3 leads to a corresponding 
change in the exponents Ai for the parameters. In this case we 
obtain Ai = 1.8 ± 0.1, A2 = 0.24 ± 0.02, A3 = 0.23 ± 0.02, 



A4 = —0.77 ± 0.05. CoiTesponding results were found for 
several other exponents for the critical shear, see Fig. |6] This 
clearly illustrates how the dependence of the steady state on 
the particle and flow properties is influenced by the frag- 
mentation mechanism. We find that a fragmentation rule of 
Scrit — ^01^^ leads to a relationship of the form of Eq. 
where the scaling exponents are linear functions of j. 

We also mention that in addition to the steady state dis- 
cussed here the transient behavior of the system is greatly 
influenced by the particle and flow properties. One of the rel- 
evant quantities for this transient is the time it takes for the 
system to reach the steady state, which is for example very 
important in the formation of rain in clouds. This time de- 
creases strongly with an increase of the coagulation rate, for 
example due to an increase of shear or collision efficiency. 
An increase in fragmentation, for example due to increased 
fluid shear or coagulate stren gth, also decreases the time to 
reach the steady state (see also lZahnow et al. How- 
ever, a detailed discussion of these transient times is beyond 
the scope of this paper and will be published elsewhere. 

3.3 Size distributions 

Finally, we discuss the size distributions of the coagulates 
and illustrate that their shape depends on the fragmentation 
process. Let us first look at the size distribution in the case 
of size-limiting fragmentation for the same parameter values 
as in the previous part. 

In this case the size distribution is fairly broad, covering 
almost the complete range between the smallest size and the 
maximum allowed size distribution (dashed line in Fig. Ilfa)). 
Comparison with other simulations for different fragmenta- 
tion mechanisms shows that the shape of the size distribu- 
tion depends greatly both on the number and on the size 
distribution of the fragments that are created during frag- 
mentation. This was already noted for fractal aggregates 
in Zahnow et al. ( 2009a) where varying size distributions of 
fragments, for example fragments of similar size or very dif- 
ferent size, were studied in the case of binary fragmentation. 
Here, we also compare our results with binary, ternary and 
quarternary fragmentation, i.e. the creation of two, three or 
four fragments instead of the lognormal distribution of the 
number of fragments described in the previous section (see 
Fig. [8]i. In particular, the width of the distribution is greatly 
influenced by this change in the number of fragments. Binary 
fragmentation leads to a single, sharp peak. Ternary frag- 
mentation leads to two broader peaks in the distribution and 
for quarternary fragmentation three peaks can be seen. These 
peaks merge into a broad plateau if the number of fragments 
is not deterministic but instead can vary as is the case for the 
lognormal fragmentation mechanism. 

In the case of shear fragmentation the situation is slightly 
different. The distribution is fairly broad, but with a long ex- 
ponential tail towards larger size classes (Fig. EJb)). The fig- 
ure clearly shows the increase in the width of the distribution 
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Fig. 7. Size distributions of the number of primary particles in co- 
agulates, normalized by the total number of coagulates: (a) size- 
limiting fragmentation for two different values of the maximum 
stable coagulate size amax and (b) shear fragmentation, for two 
different values of the coagulate strength. 
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Fig. 8. Size distributions of the number of primary particles in 
coagulates, normalized by the total number of coagulates for size- 
limiting fragmentation for a critical size of amax = 60 for different 
distributions of the number of fragments, two (binary, solid line), 
three (ternary, dash-dotted line) and four (quarternary, dashed line). 



with increasing coagulate strength, that was also indicated 
by the standard deviation (compare Fig. However, here 
it becomes obvious that the statistics degrades rapidly with 
increasing 7, as the number of coagulates available in the 
system is decreasing and hence the exponential tail becomes 
less visible. Especially in the tails of the distributions fluctu- 
ations become very large, as only a few coagulates of these 
sizes exist at all. For shear fragmentation different numbers 
of fragments, for example binary, ternary or quarternary frag- 
mentation do not influence the shape of the size distribution, 
only the mean of the size distribution is shifted towards lower 
values for an increasing number of fragments. 

The exponential tail of the size distribution is a feature that 
has also been observed for coagulatio n and fragmentatio n 



of marine aggregates in tidal flats (e.g. iLunau et al. ■ I2006h . 
A num erical comparison of d ifferent fragmentation mech- 
anisms jZahnow et al. ■ l2009ah showed that this exponential 
tail is a typical feature, when the coagulates are assumed to 
break into fragments of very similar sizes. However, differ- 
ent distributions of the fragments, for example an erosion- 
like process where some very small and some larger frag- 
ments are created lead to different size distributions of the 
coagulates. 

We note that the shape of the size distributions remains 
constant when the particle and flow properties, e.g. the co- 
agulate binding strength 7 or t he collision efficiency -Xc are 
varied. As already shown in IZahnow et al. I (l2008l l2009bl) 
the size distributions will collapse into each other for differ- 
ent parameter values when rescaled with the mean coagulate 
size. 

We also note that in the parameter ranges studied here, 
the tails of the size distributions consist of coagulates with 
Stokes numbers of order 1 (a ^ 100), which is of the order 
of the correlation time of the flow. This stron gly affects the 



clustering properties of these particles (see e.g. lCencini et al 

iBec et al.L l2007h However, since in our case there are 
typically only a very few particles in this range of Stokes 
numbers, this is not expected to significantly influence the 
properties of the steady state size distributions which we fo- 
cus on. 



4 Scaling behavior of the average coagulate size 

In this section we try to give some theoretical insight into 
the scaling of the average coagulate size in the steady state. 
We emphasize that it is not possible to find a closed equation 
for this average size of coagulates in the steady state, mainly 
because expressions for the collision rates and fragmentation 
rates are not known. 

The question of the collision rates of inertial particles is 
a topic of ongoing research, and wh ile some advances have 



been made (see e.g. lBec et aliuOOSh . even for the fairly sim 



pie flow situation used in this paper no equations exist that 
could be applied. For the fragmentation rates, the situation is 
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slightly different. While the fragmentation rate correspond- 
ing to the fragmentation model used here is known in princi- 
ple, the integral involved can not be solved in closed form. 

In the following, we will make some approximations about 
the collision rates and fragmentation rates and try to derive 
an equation for the average coagulate size in the steady state. 
While this approach makes it clear how a scaling of the av- 
erage size with particle and flow properties comes about, the 
result from the particle based model can not be fully recov- 
ered. This derivation also illustrates many of the difficulties 
associated with a mean field approach for coagulation and 
fragmentation, which do not appear in the particle based ap- 
proach described in this work. 

4.1 Equation for the average size 

To estimate the scaling behavior of the average number of 
primary particles in the steady state we start by assuming 



N{a,t) 



of 



a continuous probability distribution p{a, t) — -^vXtJ 
the coagulate^ Here, N{a, t)da is the number of coagu- 
lates consisting of a number of primary particles in the range 
[a, a + da] at time t and N{t) — daN{a, t) is the total 
number of coagulates at time t. 

The average of a with respect to p is then defined as 



(a) := daap{a,t)— da a 



N{a,t) 



(8) 



From this definition, we can derive the equation for the evo- 
lution of (a) by taking the derivative with respect to t. We 
obtain 



d 
di 



iN{a,t) 

, oo 



N{t) J N{t) 





(9) 



Defining the relative rate of change /j,(a, t) as 
p,{a,t) 



iN{a,t) 



N{a,t) ' 
equation (|9]l reduces to 

d 



dt 



(a) = {p.a) - (a) {^i) 



(10) 



(11) 



Using this equation the average number of primary parti- 
cles can in principle be calculated for all times, if the rela- 
tive growth rate fi{a, t) is known. However, this quantity is 
exceedingly difficult to determine and to date no complete 
derivation of yLt(a, t) even for very simple cases has been 



' In our model the distribution is in fact discrete, but the resulting 
sums can not be evaluated analytically. 



found and only some approximations are known. We will 
see in the next part that one of the reasons for this difficulty 
is that /i(Q!, t) depends on many properties of the whole sys- 
tem, such as the full probability distribution p(a, t). 

We emphasize that this is one of the key advantages of our 
individual particle based approach, as it only requires knowl- 
edge of the properties of the individual particles and not of 
the whole system. In addition, our approach can be used 
for the numerical calculation of the relative growth rate and 
other global quantities if the individual particle properties are 
known. 

4.2 The relative growth rate 

The equation for the relative growth rate fo r coag ulation 
Mcoag(a,i) was developed bv ISmoluchowskil (119171) . For 
each value of a there is an increase in N{a^ t) due to smaller 
particles coagulating so that their combined size is a and a 
decrease in N{a,t) due to particles of size a coagulating 
with any other particles. Formally, this can be written as 



1 



^ I'da' (x(a',a-a')iV(a',t)- 

(12) 



N{a,t) 

■N{a ~ a' ,t)C{a' , a - a' ,t)^ 
da' {x{a' ,a)N{a' ,t)N{a,t)C{a\a) 



where x(q^'i ct) is a collision efficiency, i.e. the probability 
to coagulate upon collision and C(a',a,i) is the collision 
kernel, i.e. the collision rate between particles of size a and 
a' at time t. Similarly, a growth rate due to fragmentation 
/Xfrag(Q:, t) can be developed, where N{a, t) increases due to 
larger particles breaking up so that the fragments are of size 
a and N{a, t) decreases due to particles of size a breaking 
up. This leads to 



Mfrag(a,i) 



1 



N{a,t) 



da 



' (^'d{a',a)N{a,t)F{a',t) 



N{a,t)F{a,t) 



(13) 



where i?(a', a) is the probability that a coagulate of size a' 
leads to a fragment of size a when it breaks, -d contains there- 
fore the information about the number and size distributions 
of fragments, e.g. binary or ternary fragmentation. F{a,t) 
is the fragmentation kernel, i.e. the fragmentation rate for a 
particle consisting of a primary particles. 

Both the coagulation and fragmentation kernels will in 
general not only depend on a and a' but also on system para- 
meters, for example the coagulate strength or the turbulence 
level in the flow. The total relative growth rate is then given 
by //(a, t) = //coag(a, t) + Aifrag(a, t). 
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In recent years much effort has gone into finding approx- 
imations for the collision kernels. But in particular when 
particle inertia plays a role, effects such as preferential con- 
centration and the occurrence of caustics lead to drastic mod- 
ificati ons of the collisio n kernels that are still not fully under- 
stood teec et al.ll2005l) . 



The fragmentation kernel poses a very different problem. 
On the one hand it seems to be easier because it only involves 
individual coagulates. On the other hand it can be extremely 
complicated because the microscopic properties of the coag- 
ulates play a very important role and generally both the frag- 
mentation kernel and the distribution of fragments are not 
well understood. 

4.3 Estimating the scaling of the average coagulate size 

It is clear from the previous section that estimating the rel- 
ative growth rate and with that the average coagulate size 
requires some information about collision and fragmenta- 
tion kernels. We mention again that particularly for estimat- 
ing these quantities and for comparing them with theoretical 
predictions of these quantities our individual particle based 
modeling approach is most useful. 

The size distribution of the coagulates in the steady state 
(see Fig. |7]i can be well approximated by an exponential. 
Therefore we take N{a, t) = :^e-"/<">. 

As a first approximation for the coagulation, we assume 
that differences in radius between particles are small. Then 
the problem of the collision rates reduces to that of two par- 
ticles with the same Stokes number, which is given by the 
mean Stokes number of the two particles. Since there is 
no interaction between particles through the flui d the colli 



sion k ernel in our model is then approximated by dBec et al 
2005h . 



C(a, a') cx ei/2(r(a) + r(a'))*'<^*''^^ 



(14) 



where e is the average dissipation of turbulent kinetic energy 
in the flow and the exponent is a function of the mean 
Stokes number (Stri) of the two particles. r{a) is the ra- 
dius of a particle consisting of a primary particles, here this 
is given by r{a) cx a^/'^. No analytical expression is known 
for 9{{Stri)), only two limit cases. For no inertia the expo- 
nent is 3 and for (Strj) — > oo the expon ent approaches 2. 
Numerical results e.g. bv lBec et al. I (120051) suggest that 6 de- 
creases monotonically for increasing (Strj), but no explicit 
equations are given. 

To illustrate in principle the calculation of the average 
size we concentrate on the limit case of no inertia, where 
9((Str,)) = 3. This is the so-called rectilinear shear kernel 
(iThomas et allll999l) . ]n addition, we use a constant colli- 
sion efficiency, i.e. x(q:', a) = Xc- 



Using these approximations, the relative growth rate due 
to coagulation in Eq. (fTSl i can be calculated. We obtain 



Mcoag(a,i) = Ci 



9 + 47rV3 

c 

18 



(a)^-2r(|)(a)^/^aV3 

3r(|) 



(15) 



where ci is a proportionality constant. 

For fragmentation, the first approximation for the proba- 
bility i}{a' , a) is that ^{a' , a) — 2S{a' — 2a), where S{x) is 
the Dirac delta function. This is the case of binary fragmen- 
tation, where both fragments are of the same size. 

In our model all coagulates of the same size have the same 
critical shear Sc{ct, 7). Here we studied the case Sc — ja^^, 
where ^ > and in particular the case of ^ = 1 /3, see Eq. 
©. If the fluid shear S = (25'^ S'.y■)l/^ where Sij is the 
rate-of-strain tensor in the flow, at the position of the co- 
agulates exceeds this critical shear Sc they fragment. The 
probability for fragmentation of a given size is then only de- 
termined by the probability distribution of the shear p{S, t), 
the influence of individu al particle properties for th e frag- 
mentation kernel (see e.g. iRuiz and Izquierdo is not 
considered. Again neglecting inertia effects and assuming 
a homogeneous distribution of the particles in the flow, the 
fragmentation kernel, i.e. the fragmentation rate is given by 



/ dSp{S,t)/T{S,t) 

■Sc(a,7) 

/ dSp{S,t) 




(16) 



whe re t(S, t) is the characteristic time of the shear S, see 
e.g. iBabler et a l. (2008). In our case, F{a, t) can not be de- 
termined analytically. Generally, ^(5*, t) is a function of four 
(or nine, in three dimensions) random variables Sij. Even in 
the simple case of independent normally distributed random 
variables whic h we have here, p(S , t) can not be calculated. 

For larger a iBabler et al. I (l2008 l) argued that the fragmen- 
tation kernel is approximately given by a power law function. 
They estimated that the fragmentation kernel can be approx- 
imated by F{a,t) sa ^26^, with some constant £2- Sim- 
ilar power-law approximations for the f ragmentation kernel 
have b een found in other cases, see e.g. Ruiz and Izquierdol 
(Il997h . We will therefore continue our calculation using this 
expression. 

We emphasize that the relationship between F and its ar- 
guments a and 7 depends on the specific form of the stabil- 
ity condition for fragmentation, i.e. the specific model for 
Sc- It is through this dependence that the exponent of f of 
the stability condition Eq. ^ appears in the final result. For 
Sc — 70;^^ we obtain 
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Wrag(a,0 = C2e7" 



(17) 



These approximations for the relative growth rate can then 
be used in Eq. (fTTT l. To find the scaling behavior in the steady 



state we set 4; (a) = 0, which leads to 



dt 



^/^N{a) cx ej-^ (a) 



(18) 



where N = N{t)- (a) is the total number of primary particles 
in the systerr0. The terms on the left side are the contribution 
from the coagulation Eq. (flST l and the terms on the right 
hand side follow from the fragmentation Eq. (fTTI l. Solving 
Eq. ([T8]l for (a) leads to 



(a) cx 7i/«xy^^^/^e"* 



(19) 



We find an equation for the scaling of the average coagu- 
late size in the steady state as a function of the particle and 
flow properties. While the scaling for 7 and e is similar to 
what was found in our numerical simulations using the indi- 
vidual particle based model (see Eq. (|7]i and Fig. |6l), the scal- 
ing with Xc and N is not entirely correct. This can therefore 
not be explained fully with the approximations made here. 
However we do find that the scaling exponents depend on 
l/<^, as was found in the numerical simulations. 

However, the calculation in this section illustrates that the 
average coagulate size can indeed be expected to scale with 
the particle and flow parameters and also makes it clear how 
the dependency of the scaling exponents on the exponent 
of the stability condition Eq. ^ appears. Additionally, this 
calculation illustrates the special role of e that affects both the 
coagulation and the fragmentation in the system. It is rather 
remarkable that even though the analytical calculation of the 
scaling is only possible in a very simplified case, a similar, 
simple scaling of the average coagulate size can be found 
numerically for the full individual particle based model. 



5 Conclusions 

In the present study we described results from a coupled 
model for advection, coagulation and fragmentation of indi- 
vidual inertial coagulates. The model represents an approach 
to bridge the gap between the mean field theory that is com- 
monly used to describe larger coagulation and fragmenta- 
tion systems and a full simulation of a multiphase flow. Full 
hydrodynamic simulations of coagulation and fragmentation 
are computationally limited to systems with very few parti- 
cles and are therefore not appropriate to describe large-scale 
processes such as initiation of rain in a cloud. Mean field 



This follows from the assumption of an exponential size distri- 



bution. 



models on the other hand are capable of describing coagu- 
lation and fragmentation on such scales, but rely on many 
approximations and parameterizations. 

Our individual particle based approach was used to gain 
insights into the principle behavior of coagulates under dif- 
ferent fragmentation mechanisms and to study the depen- 
dence of the steady state of the coagulates on particle and 
flow properties. We used synthetic turbulence in the form of 
a smooth random flow to approximate the motion of particles 
in a turbulent flow, focusing on processes which take place 
below the Kolmogorov scale. Even though not all features 
of turbulent flows are captured, the results are expected to 
remain qualitatively similar in more realistic flows. In realis- 
tic turbulent flows clustering and collisions between particles 
may depend on non-Gaussian statistics and intermittency in 
the velocity field, as well as the Reynolds number and could 
also be affected by clustering at an inertial range, where the 
velocity field is not smooth. However, as long as the system 
is well mixed, we do not expect a strong qualitative change. 
The same is true for the extension to 3-d flows, where coag- 
ulation slows down, due to less frequent collisions, thereby 
mostly affecting the time scale of the approach to a steady 
state. 

The applicability of the model used here to more realistic 
problems is limited due to the computational restriction of 
the number of primary particles. However, it is well suited 
for small systems and principal studies of underlying mecha- 
nisms. A great advantage is that an individual particle ap- 
proach can easily incorporate experimental results and re- 
sults from full hydrodynamic simulations to calculate aver- 
age quantities such as collision or fragmentation rates which 
can then be incorporated into larger mean field models. 

In this work we numerically studied the steady state that 
results from a balance between coagulation and fragmenta- 
tion. Mainly, we examined average quantities that character- 
ize the steady state, such as the average number of primary 
particles per coagulate. We compared two different frag- 
mentation mechanisms, size-limiting fragmentation which is 
motivated by the hydrodynamic al instability of large drops 
settling under gravity and shear fragmentation, where parti- 
cles break due to hydrodynamic shear forces. For both size- 
limiting and shear fragmentation the transient behavior of the 
system is strongly influenced by the particle and flow proper- 
ties. In particular, enhanced collision rates, for example due 
to increased shear or increased collision efficiency greatly 
decrease the time it takes to reach the steady state. 

For size-limiting fragmentation this steady state shows few 
fluctuations and almost no dependence on the particle or flow 
parameters. The main parameter that determines the coagu- 
late size distribution in this case is the maximum stable coag- 
ulate size. The size distribution in this case is very broad, and 
covers almost all the available coagulate sizes. Different size 
distributions can appear if the number of fragments is cho- 
sen differently, for example if fragmentation is binary. The 
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shape of the size distribution is then related to the number of 
fragments that are created during fragmentation. 

The size distributions for shear fragmentation have a sin- 
gle peak with an exponential tail. This is a typical feature of a 
fragmentation mechanism whe re coagulates break in to sim- 
ilar sized fragments (compare Zah now et al.l ( l2009al) ). For 
shear fragmentation strong fluctuations in the average num- 
ber of primary particles per coagulate due to statistical fluc- 
tuations of the carrier flow appear In this case, both the 
average number of primary particles per coagulate and the 
standard deviation of the coagulate size distribution in the 
steady state change strongly with the particle and flow prop- 
erties. Simulations showed that the variation of each para- 
meter within a reasonable range is approximately indepen- 
dent of the values of the other parameters. For variations of 
the coagulate strength 7 the scaling relationships for both the 
average and the standard deviation can be inferred from the 
fragmentation mechanism. Scaling relationships for varia- 
tions of the volume fraction, the collision efficiency and the 
dissipation of turbulent energy in the fluid were derived from 
the simulation results. For each of these parameters we find 
a power-law dependence, where the exponents appear to be 
closely connected to the shape of the stability condition for 
fragmentation. We illustrated this by showing how an equa- 
tion for the average coagulate size in the steady state can be 
derived. From this we calculated scaling relationships for 
the average coagulate size in the steady state using severe 
approximations. This calculation also clarified how our in- 
dividual particle based approach can be connected with the 
mean field theory that is commonly used to describe larger 
coagulation and fragmentation systems. However, this ap- 
proach requires expressions for the collision and fragmen- 
tation rates as well as some knowledge of the coagulate size 
distribution. By contrast, our individual particle based model 
only requires knowledge if the individual particle properties, 
which turns out to be a great advantage of our approach. It 
can therefore be a very useful tool, both for obtaining esti- 
mates of global quantities such as collision and fragmenta- 
tion rates and as a comparison for results from mean field 
models. 

Our results emphasize the great importance of the frag- 
mentation mechanism for the final size distribution of coagu- 
lates in the steady state. As a consequence it is very desirable 
to design experiments to investigate the fragmentation of par- 
ticles in different applications. 

In general, the dependence of the average quantities as 
well as the size distributions on the particle and flow prop- 
erties can change quantitatively for different fragmentation 
mechanisms, in particular for different number and size dis- 
tributions of fragments created during fragmentation. How- 
ever, the qualitative picture that has emerged can be expected 
to remain the same. 
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